ELLIPTIC CURVES AND A NEW CONSTRUCTION OF 
INTEGRABLE SYSTEMS 



VLADIMIR DRAGOVIC AND BORISLAV GAJIC 

Abstract. A class of elliptic curves with associated Lax matrices is 
considered. A family of dynamical systems on e(3) parametrized by 
polynomial a with above Lax matrices are constructed. Five cases from 
the family are selected by the condition of preserving the standard mea- 
sure. Three of them are Hamiltonian. It is proved that two other cases 
are not Hamiltonian in the standard Poisson structure on e(3). Integra- 
bility of all five cases is proven. Integration procedures are performed 
in all five cases. Separation of variables in Sklyanin sense is also given. 
A connection with Hess-Appel'rot system is established. A sort of sep- 
aration of variables is suggested for the Hess-Appel'rot system. 
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1. Introduction 

Starting from a class of elliptic curves we construct a class of systems on 
the Lie algebra e(3) by the Lax representation. Equations of motion depend 
on an arbitrary polynomial a and have the form: 

x = {x, Hi} + a{x, H2} 

From the condition that the standard measure is preserved, we obtain five 
choices for a. In three cases, the systems are Hamiltonian. However, in 
two other cases they are not Hamiltonian in the standard Poisson structure 
on e(3). Integrability in all five cases is performed. For all five cases clas- 
sical integration procedure is presented. For the first three cases we gave 
an algebro-geometric integration procedure also. Using the Lax matrix, a 
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separation of variables in Sklyanin sense in the first three cases is shown. 
At the and, a connection with the Hess-Appel'rot system is established. A 
sort of separation of variables is suggested. 

2. Lax representations and elliptic curves 



Let us start from a class of elliptic curves 

2 a b c 

H =-p--- — - — 

They have the following representation: 



(1) 



d 

A 4 "' 



(2) 

where 



W 2 (A) A(A)A*(A) 



A 4 



A 1 



w(A) = -i(qX 2 +x 1 X + y 1 ), A(A) = xX + y, A*(X)=xX + y. 
Such curves appear as spectral curves of Lax matrices of the form 



(3) 



L(X) 



A* (A) 



We will consider that corresponds to real coefficients a, b, c, d. 

Let the functions M%, M 2 , M3, Fx, T 2 , T3 be generators of Lie algebra e(3) 
with Poisson structure defined by the relations 

(4) {M i ,M j } = -e ijk M k , {M i ,T J } = -e ijk T k , i,j,fc = 1,2,3. 

Assume the following change of variables: 



V 

y 
yi 



i=(/3M!-aM 3 



— /=(/3ri - aT 3 - iT 2 ), a 

-p(/3Ti - aT 3 + iT 2 ), x = -=(pMi - aM 3 + iM 2 ), 
v 2 v 2 



aTi + f3T 3 , x x = olM x + (3M 3 , q = I 2 



Xq+ Zq, 



where a 



and xq, zq, I 2 are constants. 



V x o+ z o '_ _ V x o+ z o ' 
In terms of x,y,x,y,xi,yi the Poisson structure (H]) has the form 

{x, y} = 0, {x, y} = 0, {x, xx} = ix, {x, x x } = —ix, {y, x t } = iy, 

{y, xx} = -iy, {yi,xx} = 0, {x, y x } = -iy, {x, yx} = iy, {yi,y} = 0, 

{yi, y} = 0, {x, x} = -ixx, {y, y} = 0, {x, y} = -iyx, {x, y} = iyx- 

To each elliptic curve ([2]) we correspond a family of dynamical systems: 



(5) 



L(A) 



1 

2h 



L(X), 



X 2 L(X) - a 2 L(a) 



X 



The matrix L(X) in ([5|) is of the form (|3|) and a is an arbitrary polynomial 
in generators of algebra e(3). 
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Observe that matrices L given by © satisfy 



(6) 

where 



L(\),L(y) 



L{\) = L{\) ® 



r(\-»),L(\)+L(») 



1 
1 



, L(jj) 



1 
1 



with permutation matrix as an r-matrix 

"1 0" 



r<A) = ± 



10 
10 
1 



3. Equations of motion 

The aim of this section is to provide the initial analysis of dynamical 
systems ([5]). 

The Poisson bracket (j4]) is degenerative. As it is well known, there are 
two Casimir functions: 



(7) 



F x = MxTx + M 2 r 2 + M 3 r 3 , F 2 



r? + rl + r§. 



Thus, a symplectic leaf, defined by conditions Fx = cx,F 2 = c 2 is a four- 
dimensional manifold. For integrability in Liouville sense on e(3), one first 
integral more beside the Hamiltonian is necessary. On the other hand, if a 
system is not Hamiltonian, generally speaking, five first integrals of motion 
for integrability in quadratures are required. But, if a nonhamiltonian sys- 
tem has an invariant measure, then, according to the Jacobi theorem, for 
integrability one needs only four first integrals of motion. 
Equations ([5]) can be rewritten in the form 



z T 2 + az M 2 , 

x T 3 - z T 1 + a(x M 3 

-x T 2 - ax M 2 , 

+ a2: r 2 , 
+ a(x T 3 

- ox r 2 . 



ZqT 



Mi 

M 2 = x r 3 - z T 1 + a(x M 3 - z Mx), 
M 3 

(8) fi 

r 2 
r 3 

We have the following Proposition. 
Proposition 1. System ([8]) can be rewritten as: 
M i = {M i ,H 1 } + a{M i ,H 2 }, 
t i = {T i ,H 1 } + a{T i ,H 2 }, i 



r 2 M 3 




r 3 M 2 




h 




r 3 Mi 




riM 3 




h 




r!M 2 




r 2 Mi 


h 



(9) 



1,2,3, 
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where 

#1 = — 7TT + («ori + Z r 3 ), #2 = X Mi + Z M 3 . 

^2 

For a general polynomial a, the system ([8]) is neither Hamiltonian in the 
Poisson structure (J4j), nor preserves the standard measure. Simple criterium 
for preserving the standard measure is given by: 

Proposition 2. The system ([8]) preserves the standard measure if and only 
if the polynomial a satisfies the condition: 

{a,x M 1 + z M 3 } = 0. 

Proof. The divergence of the vector field in (jSJ) is equal to {a, xqM\ -\-zqM 3 }. 

□ 

As a consequence of Proposition [2] we have: 

Proposition 3. In the following five cases, the standard measure is pre- 
served 

(i) : if the polynomial a is a Casimir function: a = M\Tx + M 2 T 2 + 
M 3 r 3 ; 

(ii) : if the polynomial a is a Casimir function: a = rf + r| + r§; 

(iii) : if a = xqM x + z M 3 ; 

(iv) : ifa = x T 1 + z T 3 ; 

V 



(v): ifa = Mf + M 2 2 + Ml 



Theorem 1. If xq ^ 0, or zq ^ 0, in the first three cases given above, the 
systems are Hamiltonian, while in the fourth and the fifth case, the systems 
are not Hamiltonian in the Poisson structure 

Proof. If a is a Casimir function, for an arbitrary function / we have 

{/, Hi} + a{f, H 2 } = {/, H x + aH 2 }. 

Hence, in the first two cases the systems are Hamiltonian with Hamiltonian 
functions 

H = H l +aH 2 . 
In the third case, since a = H 2 , we have 

{x\ #1} + H 2 {x\ H 2 } = {x\ H x + ^}, 

where x l , i = 1,...,6 are coordinates M\, M 2 , M 3 , Y\, T 2 , T 3 . Thus, the 
system is also Hamiltonian with the Hamiltonian function 

The fourth and the fifth cases are more complicated. For a system (|8|) to be 
Hamiltonian, the polynomial a needs to satisfy 

a{x\H 2 } = {x\H} 
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for some function H. We get the condition 
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(10) 



^ i ( dH 2 dH\ 
> \x ,x J \ \ a——. — — -r = 0. 

7 = 1 



In the fourth case, when a = xqTi + zoT 3 , the system of linear partial 
differential equation (fT0|) becomes 

dH dH dH dH 

M3 dW 2 ~ M2 m 3 + XoZoTlM2 + z2 ° T3M2 + r3 w 2 " T >m; = °' 

dH dH 

Ml dM 3 ~ ^dMl + xlTlM * + X0Z ° T3M3 

dH dH 

- xozqTxMx - z 2 r 3 M! - r 3 — + ri— = o, 

Oil OL 3 

dH dH dH dH 



dH „ dH 



T3 dW 2 - T2 dW 3 + X0Z ° rir2 + = °' 



dH „ dH 



r i7TiT~ ~ x o r i r 2 - x z T 3 T 2 = 0. 



dM x 8M- 

It is easy to see that the system (fTTj) consists of only four independent 
equations 



(12) 



dH 

2 dM l 


= v x 


dH 

dM 2 


■ + xgrir 2 + x z r 2 r 3 , 


r dH 

2 dM 3 


= r 3 


dH 

dM 2 


■ + xo^rii^ 


+ ^r 2 r 3 , 


^dT, 




dH 
dT 2 




r i 5 # 

r 2 2 dM 2 


dH 

T2 dv 3 


= r 3 


dH 
dT 2 




T3 M 2 dH 

r 2 2 dM 2 



From the first and the third equation (|12p . one gets 

d 2 H d 2 H _ 2 

aridMi ~ aMiari ~ x °' 



If xq 7^ 0, the condition ([10]) is not satisfied. Consequently, the system is not 
Hamiltonian. In the case when zq ^ one gets contradiction by observing 
that 92 H - 92 H ^ n 
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Following the same procedure given below, in the fifth case, when a 
Ml + M 2 2 + Mf system ([TO} becomes 
(13) 

M 3 ^- - M 2 ^- + z M 2 (M 2 + Mf + Mf ) + T 3 |^ - T 2 |^ = 0, 
A dM 2 dM 3 u K 1 2 3J 6 dT 2 dT 3 

dH OH 

Ml dW 3 ~ Mz dM[ + X ° M ^ M " + M * + M ') 

flu fl rx 

- zqM^mI + M| + Mf) - r 3 — + Tx— = 0, 

M 2 ^- -M x ^-- x M 2 (M 2 + Ml + Mf ) + T 2 |^ - = 0, 

dM x dM 2 u A 1 2 37 sri dr 2 

flH BH 

T3 dM 2 - T *im + z ^ iM " + Ml + M * ] = °' 

d H BH 

Tit^ T 3 - + (x r 3 - z T 1 )(M 2 + M? + M 3 2 ) = 0, 

BM 3 6 dM x V U 3 U A 1 2 37 

^ " " X ° r2(M ' + M2 + = 

System (JT3J) has four independent equations 



(14) 



r dH - 

2 BM 1 " 




dH 

0M 2 


■ + x r 2 (M 2 


+ M 2 2 + Mf) 


d# 

2 <9M 3 " 


r 3 


dH 

0M 2 


■ + z r 2 (M 2 


+ Mf + Mf ) 






dH 
dT 2 


+ Mi 

9M 2 


r 2 2 9M 2 ' 


r2 sr 3 " 


r 3 


dH 
dT 2 


+ M 3 - 

6 dM 2 


r 3, r d# 

r 2 2 9M 2 



From the first and the third equation ()14|) . one get 
d 2 H d 2 H M 



2x ^(M 1 r 2 -M 2 r 1 ). 



dTxdMi dMidT! T 2 
Consequently, in the fifth case system is not Hamiltonian. □ 

Regarding integrability of given five cases, we have simple Proposition. 

Proposition 4. 

(a) A function F is a first integral of equations ([8]) if it satisfies 

F = {F,H 1 } + a{F,H 2 } = 0. 

(b) The Casimir functions F\ and F 2 and functions H\ and H 2 are integrals 
of system flSJ) for any polynomial a. 

Finally we have: 

Theorem 2. The system ([8]) in cases (i)-(iii) is completely integrable in 
Liouville sense. In cases (iv) and (v), system (J8|) is integrable in quadratures. 
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4. ALGEBRO-GEOMETRIC INTEGRATION PROCEDURE OF THE SYSTEMS 

We pass to an algebro-geometric integration procedure, based on con- 
struction of the Baker- Akhiezer vector-function. For more details about the 
Baker- Akhiezer functions see 0[8j|9]. 

The dynamical systems have Lax representation 

L(X) = [L(X),A(X)], 

with L(X) given by © and 

A 2 L(A) - a 2 L(a) 



A(X) 



2/ 2 (A 



The corresponding spectral curve T is elliptic curve ([2]). 
As usual, we consider the following eigenvalue problem 

(15) (^«)*< ( > P > = °- 

L(A)*(t, P) = ^(t, P). 

with a normalization 

(16) ^(O,^) + ^ 2 (0,P) = 1, 

where P = (A, fj,) is a point on the spectral curve T. 

Let us denote by $(t, A) the fundamental solution of the system 

! + W)*( t ,A) = o, 

normalized by the condition <£(0, A) = 1. 

Let us denote by 00 + and 00 ~ the two points on the curve T over A = 00, 
with /i = H2\/xq + Zq and [i = — H2\J x\ + z 2 respectively. 

Proposition 5. // the polynomial a is a first integral of motion, then the 
vector- function ^(t,P) satisfies the following conditions: 

(a) In the affine part of the curve T, the vector-function ^>(t, P) has two 
time independent poles, and each of the components 1 & 1 (t,P) and 
^ 2 (t,P) has one zero. 

(b) At the points 00 + and oo~ , the functions X I /1 and ^ 2 have essential 
singularities with the following asymptotics: 



y 2 {t,p) 



i(v^I+^ (A+a)+f )t (pi 1 )) ^ 



P^OQ- 
P^OO + 

p OCT 
p OCT 
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(c) The asymptotics have the form 



&(t,P) = e ~^ x ° +z ° (A+a)+ t)* 



1 



/ 2 vV4+^o A 



I 2 V2y/x% + z 2 Q A 



+ 0(1/A 2 



+ 0(1/A 2 



P 



Proof. Since 

*(t,P) = *(t,A)*(0,P), *(0,A) = 1, 

the poles of P) coincides with poles of ^(0, P) (<3?(i, A) is holomorphic) , 
so they are time-independent. Using normalization (|16|) . and (|15|) one can 
calculate that \&(0, P) has two poles Pi and P2 in the affine part of T. One 
can also conclude, using din and din ^2 that each of the components ^1 
and VP 2 has a one zero. This proves part (a). 

In order to find asymptotic of ^f 1 at 00 + and oo~ one needs to consider 



(mOf 1 ))' 



£1 



-A n - A 



12 



* 2 
*7 



Since |£ = ^ L ^ n = „^x 22 , and using asymptotics 

A* = i h\fxj+z$ + ^ + 0(1/A) , P^oo+, 



A 



and 

we have: 
(m^ 1 ))- 
and 



A* : 
~2 



vMg + ^o 2 + y + O(VA) 

2^2 



P -> DO~ 



s o + z o (A-a) 



, .'f, • : 2 - • 1 A). P • x 



(ln^ 1 ))' 



^x 2 + z 2 (A + a) + g+O(l/A), P 



00 



Since 



x 
x 



i<i\/xl + zl+iyj xl + zl |. 



we finally get: 
(ln^ 1 ))- = 

and 

(ln(^ 1 ))' 



^x 2 + z 2 (A + a)-g + | + 0(l/A), P 



00 



2V * 2 + z 2 (A + a) + || + 0(l/A), P 



Similarly for the function VP 2 we have: 
(ln(* 2 )y 



--\ x 



IX\ 



+ ^ 2 (A + a)-^ + O(l/A), P 



00 



00 
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and 



(ln(* 2 )y = y x l + zl (A + a) + g + ! + 0(1/A), P 
This proves parts (b) and (c). 



oo 



□ 



Remark 1. The vector- function ^f(t, P) satisfies standard conditions of the 
2-point Baker- Akhiezer function. From the parts (a) and (b) of Proposition 
[5]we can construct 9(t, P). Using the part (c), one can reconstruct x(t) 
(and x(t)) from the Baker-Akhiezer function. 

Remark 2. All formulae from the proof of Proposition [5] are still valid even 
if the polynomial a is not an integral of motion. But in that case statements 
(b) and (c) of the Proposition are not valid any more. 

Now we will give explicit formulae for the Baker-Akhiezer function in 
terms of the Jacobi theta-function 6h{z\t) with characteristics [i, i]. Ex- 
plicit formulas are similar to the Lagrange case (see 

Let us fix the canonical basis of cycles A and B on V {A ■ B = 1), and let 
uj be the holomorphic differential normalized by the conditions 



LO 



2z7T, 



LO 



T. 



I A JB 
Theta-function 9\\{z\t) is defined by the relation 



^exp 



\ T ( n + \) + (* + **■)(« + -) 



and satisfies 

U (O) = O, &ii{z + 2m) 



hi(z + r) = -exp(-z- -t)6u(z). 



Let Q + and f2 be differentials of the second kind with principal parts 

~~ I V x o "+ z o and +|\/ x o "+" z o ^ a ^ °° + an d a ^ °°~ respectively, normal- 
ized by the condition that A-periods are zero. Let us introduce differential 
0, = S7 + + . We will denote by U the P-period of differential f2, and by 
c + and c~ the constants: 

" p i I 

fl = --Jx 2 + z% \ + c + + 0(l/\), P^oo+ 

'Po 2 v 

»P „• , 

^ =- + t \ / - < 7, -!- : n A + < ■ + </( l/A). / ' — x 

'Po 



+2^0 + ^ + ^ + 0(1/^, P 



Proposition 6. JTte Baker-Akhiezer functions are given by 

6 n (A(P + oo^ 



^(tjP) = ci exp 
* 2 (t,P) = c 2 exp 



p 
Po 

(f P Q 
JPo 



I III. 



Pi - P 2 ) + tU) 



2 V 



n (A(oo+ + oo- - Pi - P 2 ) + tU) ' 

9 U (A(P + oo- - Pi - P 2 ) + tip 
5n (^(oo+ + oo- - Pi - Pa) + ttf) ' 



10 



VLADIMIR DRAGOVIC AND BORISLAV GAJIC 



where constants c\ and C2 are 

= e u {A{P - oo+))M*4(oo- - Pi))0iM(<x>~ - p 2)) 
Cl en(^(oo--oo+))e 1 i(X(P-Pi))eii(A(P-ft)) ' 
= e u (A(P - oo-))M.4(oo+ - Pi))gn(.4(oo + - P 2 )) 

02 e n (A(oo+ - oo-))e u (A(p - Pi))e n (A(p - p 2 )) ' 

and A is the Abel map, and Pi and P2 are the poles of the function ^. 
Proof. The proof is based on general theory (see E] ) • D 



5. Classical integration of the systems 



We are going to show that in all considered five cases integration of the 
systems can be reduced to elliptic integrals. 

Let us change coordinates by the rotation of xOz plane: 



Xi = aMi + (3M 3 X 2 = M 2 X 3 
Y l = aT l + [3T 3 Y 2 = T 2 Y 3 

Differential equations ([8]) of motion become 

X x = 

X 2 = y/xl + z* (Y 3 + aX 3 ) 
X 3 = - 



-/3Mi + aM 3 
-pr x + aT 3 



(17) 



Y-i 



1 

h 
1 



xi + (Y 2 + aX 2 ) 



(X 3 Y 2 - X 2 Y 3 ) 



Y 2 = f(X x Y 3 - X 3 Y X ) + a^xl + z 2 Y 3 
1 



Y 3 = jVC 2 Yy - X X Y 2 ) - ay/ 4 + zl Y 2 . 



The first integrals are 
Pi 
P 2 

(18) 

H 2 

From (|17p we have: 



X X Y X + X 2 Y 2 + X 3 Y 3 = ci 
Y* + Y 2 2 + Y£ = c 2 



xl + xl + x$ 

2I 2 y/x% + Z$ 

X x = d 2 . 



+ Y 1 = d 1 



xl + xl 



x 2 x 3 — X 3 X 2 



(xl + zl) 



x + z 



(Y 2 2 + Yi) + a\X 2 2 + Xl) + 2a(X 2 Y 2 + X 3 Y 3 ) 
X 2 Y 2 + X 3 Y 3 + a{Xl + Xl) . 
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Using the first integrals (fT8|) . one gets 

X\ + X\ = (xg + zg) [c 2 - if + a 2 (X 2 + X 2 ) + 2a( Cl - d 2 Y 1 ] 

X 2 X 3 - X 3 X 2 = ^x 2 + zl [d - d 2 y x + a(X| + X 2 ) 

Introducing polar coordinates 

X 2 = p cos o", X 3 = p sin ex 
the last equations become 



■2 



(19) 



+ p 2 a 2 = {x 2 + z 2 ) [c2 - if + a y + 2a(ci - d 2 *i)] 



;§ + z 2 [d - d 2 Y 1 + ap 2 



Let us subtract the square of the second equation multiplied by 4 from the 
first equation multiplied by 4p 2 . After simplifying one gets 
2 

= 4(x 2 + z 2 ) [p 2 (c 2 - if) - (d - ^yo 2 ] . 



Finally, using 
(20) 

and denoting p 2 

(21) 

where 



d 2 + p 2 



«, we get 



2/2^0^ 



5m 2 - Cm - D 



-4A^2 + 



1 ^2 



C = 4(x 2 +z 2 )U 2 -c 2 + 



rf 2 (ci - cfeyl) 
hy/xfTzfi 



D = 4(xq + Zq)(ci - d 2 A) z 
A = d x 



4 



2I 2 ^x 2 + z{ 
So, the following proposition is proved: 

Proposition 7. The function u(t) is an elliptic function of time. 

Let us remark that u (and consequently p) does not depend of a choice 
of the polynomial a. 

Having u(t) as a known function of time, one can find p(t) as a known 
function of time. In order to reconstruct X 2 and X 3 , one needs to find a as 
a function of time. 



12 
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From the second equation of (fl9|) . using ([20]) . we have 




) 



+ ap 2 {t) 



The right hand side of the last equation is a function of time and of the 
polynomial a. When a is a first integral of motion (in the first three cases 
defined above) e.q. when a = const, then right hand side of the last equation 
is a known function of time and one can find a by quadratures. In the fourth 
case 



So, in this case a(t) is a known function of time and one can find a by solving 
a differential equation. Similarly, in the fifth case 



is again a known function of time and a differential equation for determining 
a can be solved. Knowing p and a as functions of time, one can reconstruct 
X2 and X3. From (|20p. one finds Y\ as a function of time. Finally, using the 
differential equation for Y\ from (|17p . and the second first integral F2 = C2 
from (|18p one can reconstruct Y% and Y3. 

Two elliptic curves appeared here. The first one T, has been defined by 
the equation ([1]), and it was the curve from witch we started. The other one 
r', given by 



1 2 

corresponds to the solution of the differential equation (|2ip . A natural ques- 
tion is how these two curves are related. We have the following proposition: 

Proposition 8. The elliptic curves T, defined by the equation (pQ) and V 
defined by (]22p are isomorphic. 

Proof. A direct calculation gives that j-invariant of both curves are the 
same, so they are isomorphic. □ 



One of the oldest methods in the theory of integrable dynamical systems 
is separation of variables. Originally, this method was built in order to 
find exact solutions of the Hamilton- Jacobi equations. In the middle of 
1990's, Sklyanin in his celebrated paper [IB] , introduced a new concept of 
separation of variables which was related to modern techniques in the theory 
of integrable systems such as the inverse scattering method and the Lax 
representation. 




a = Ml + Ml + Mf = X\ + Xl + X$ = d 2 2 + p 2 (t) 



(22) 



v 2 = -^-Bu 2 -Cu-D 



6. Separation of variables 
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If a completely integrable system on 2n-dimensional symplectic manifold 
is given with n functionally independent commuting first integrals of mo- 
tion F±,...,F n , then variables Ai, A n , u\, \x n are separated if they are 
canonical e.q. 

{Ai,/Uj} = Sij, {X i: \j} = {fM,Hj} = 0, i,j = l...n, 

and if there exist n relations <E>j such that 

(23) $ i (X i , f M,F 1 ,...,F n )=0, i = l,...,n. 

The Sklyanin magic recipe gives separation variables Aj as poles of the 
properly normalized Baker-Akhiezer function. The canonically conjugated 
variables are corresponding eigenvalues of the Lax matrix L(Aj) (for details 
see [IS]). 

In the case of algebra e(3), since symplectic leaves are four-dimensional, 
one needs to find four separation variables Ai, A2, jUi 5 fJ-2- 

Proposition 9. For the Hamiltonian systems (i), (ii), (Hi) defined in Propo- 
sition^ the separation variables are 
(24) 

y f nr~i (aM!+/3M 3 )x (oTj + pT^x 2 ' 
Ai = -, / i 1 = t^xg + zg + ^ 

A 2 = x, u 2 = + /3M 3 ). 

x 

Corresponding separation relations are: 

2 w(-Ai) 2 -2A(-Ai)A*(-Ai) , 

(25) nl = -i ^ A — ^ i ^, A 2/ u 2 = -iH 2 = const. 

Proof. Following jTB| [T5] , the separation variables for the systems ([3J) satisfy 
the equation 

(26) (10)adj(L(A)-/i-l) = 0, 

which corresponds to the standard normalization d?o = (1 0). From (|26|) one 
gets only one pair of separation variables Ai,/Ui. The other two variables 
A 2 ,//2 are calculated from the asymptotics of the matrix L(X) when A goes 
to infinity. 

By direct calculations we see that variables are canonical. □ 

In terms of Hi , H2 , Fi , F2 , the separation relations can be rewritten in 
the form: 



1 



Mi = T4 ( - + zl)\\+2I 2 Jxl + zl H 2 X\ - 2J 2 #iA? + 2F 1 \ l - F 2 ), 

A i 

H 2 = 1X2^2- 
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7. Connection with Hess-Appel'rot system 

In this section we find a sort of separation of variables for the Hess- 
Appel'rot case of motion of heavy rigid body about fixed point. 
The equations of motion are the Euler-Poisson equations: 

(27) M = MxQ + Tx X , r = TxO. 

Here M is the kinetic momentum vector, £1 is the vector of angular velocity, 
r is the unit vertical vector, and x = ( s 0j 2/0 j -^o) is the radius vector of the 
center of masses. Connection between M and f2 is given by M = IQ, where 
/ is the inertia operator, and we can assume that I = diag(Ii, I2,I 3 ). 

Equations (|27p are Hamiltonian on the Lie algebra e(3) in the standard 
Poisson structure given by @ with the Hamiltonian function: 

(28) H RB = ^(M,n) + (r,n). 

Hence, for integrability in the Liouville sense of equations (I27p . one needs 
one first integral of motion more. 

The Hess-Appel'rot case of rigid body is introduced by Hess and Appel'rot 
(see [UEE]). It is defined by conditions 

Vo = 0, 

(29) 



X{) 



y/h{I 2 - h) + ZoVWl ~ h) = 0. 



Instead of the fourth integral, in the Hess-Appel'rot case there is an invariant 
relation 

(30) x M 1 + z M 3 = 0. 

In [1] (see also [5]) a Lax representation for the Hess-Appel'rot system was 
found. The spectral curve is reducible. One component is a rational curve, 
and the other one is an elliptic curve. In order to find a sort of separation of 
variables for this specific situation, we give first another Lax representation. 
It is based on following observation: 

Proposition 10. On hypersurface (I30p the equations of the Hess-Appel'rot 
system are equivalent to the Lax representation (|3|) with a = ^=£=J . 

This Lax representation for the Hess-Appel'rot system was given in [10J. 
The Lax representation of the same form for the Lagrange top was con- 
structed in [15] . 

As a consequence of Propositions [10] and [9l we have: 

Proposition 11. The separation variables for the Hess-Appel'rot system 

are 

(31) 

. y .( rY~l (aM l +f3M^)x (ar 1 +/3r 3 )x 2 

Ai = -, fii = i [I 2 \ Xq + zZ 1 ^ 

x \ v y y z 

A 2 = x, ii 2 = ~-(aM 1 + f3M 3 ). 
x 
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Corresponding separation relations are: 

(32) ^ ^f-yW , „ 2 = . 

In terms of Hi,H2,Fi, F2 separation relations can be rewritten in the 
form: 

Mi = ^4 ( - + zl)\\+2I 2 ^xl + zl H 2 \\ - 212^X1 + 2F 1 X 1 - F 2 ), 

i\2^2 = 0. 
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